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RESEARCH ARTICLES 

New Method for Characterizing 
Dissolution Properties of Drug Powders 

PETER VENG PEDERSEN 



Abstract □ Various multiparticulate dissolution models that assume 
a log-normal particle-size distribution are fitted by nonlinear least- 
squares regression to data from the dissolution of micronized glyburide. 
Estimates of parameters describing the effective initial particle-size 
distribution are obtained, together with estimates of the specific disso- 
lution rate parameter. A dissolution equation based on an ide&l, un- 
truncated, log-normal distribution, with the single particles dissolving 
according to the Hixson-Crowell cube root law, is the best model. The 
dissolution behavior of glyburide can be well described by this model in 
terms of the specific dissolution rate parameter and one other parameter 
accounting for the distribution effect The estimation of these two pa- 
rameters represents a more exact way of describing the dissolution 
characteristics of drug powders than previous approaches. The method 
should be of interest in the quality control of drugs that may cause bio- 
availability problems because of dissolution rate-limited absorption. 

Keyphrases □ Dissolution— micronized glyburide, various multipar- 
ticulate models fit to data, rate and size distribution parameters evaluated 
□ Powders — micronized glyburide, various multiparticulate models fit 
to dissolution data, rate and size distribution parameters evaluated O 
Glyburide, micronized — dissolution, various multiparticulate models 
fit to data, rate and size distribution parameters evaluated 

The dissolution behavior of a multiparticulate system 
is determined by the single-particle dissolution behavior 
and the particle-size distribution (1, 2). Dissolution profiles 
have been related to the initial particle-size distribution 
(3, 4). However, the mathematical models used to account 
for the size distribution effect were not exact but were 
based on empirical approximations. 

A previous paper (5) showed good agreement between 
the theoretical dissolution profile based on a rigorous 
mathematical treatment and dissolution data for tolbu- 
tamide. The present paper evaluates some kinetic models 
for multiparticulate dissolution by fitting them to data 
from the dissolution of micronized glyburide. This method 
represents a new and better approach for characterizing 
the dissolution behavior of powders by evaluating both rate 



and size distribution parameters The method should be 
of value in quality control measurements of slightly soluble 
pure dmj> that can c&use bioavailability problems. 

THEORY 

Consider a multiparticulate dissolution where the single particles 
dissolve according to the Hizson-Crowell (6) cube root law or the Nie- 
bergail ei al (7) square root law, which can be written in a common form 
as: 

w Um -m Um k m t (Eq.l) 

where m - 3 for the cube root law and m ^ 2 for the ^juare root law 1 . 

Let the initial effective particle diameter* oa, be log normally dis- 
tributed, i.e., have a frequency distribution that »a well approximated 
by t he following density distribution (2): 

(Bq .2> 

•/tea — 



for j* - 1, < hi oo £ M + jo where i and j aie the lower and upper trun- 
cation parameters (Fig. 1), and Wis the normal density function defined 

by: 

N(w)~^e-i(^y (Eq.3) 

The dissolution of such a multiparticulate powder is described by 

(2): 
W 
Wo 



£ ( m ) (-K m i )<»-«> 
\ at m I \ a m f 



X — 



(Bq.4) 



< The positive constants a 3 (m - 3) and k t (m - 2) are used for s 
P W the coefficients of time ihal contain quantities such as density, t 
^oefnaeni, tht «hapt fv'uw. which are considered constant 
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and: 



V-io U 

In (PARTICLE DIAMETER) 
Figure I— Illustration of the initial (t = 0) particle-size distribution 
parameters, n, a, i, and), of a lag-normal powder. The approximating 
distribution function shown is defined in Eq. 2. The initial diameters, 
ao, are diameters of hypothetical spherical particles that approximate 
the dissolution behavior of the real nonspherical particles. The di- 
mensionless distribution shape parameters, i, j, and a, can be estimated 
from the dissolution profile (Figs. 3-6) by nonlinear least-squares re- 
gression analysis (Table I). The scale parameter, u, cannot be deter- 
mined. 

where W/W 0 is the fraction of undissolved powder at time t , F is the 
cumulative standard normal distribution function defined by: 

F(x) = -Lz C e-***dx (Eq.5) 

Ti and T 2 are given by: 

7, = n - ia for j In (K m t) < » - (Eq. 6) 

Fi-^ln (K m t) (otj\n(K m t)^u-ic (Eq.7) 

T 2 =n+ja for^ln(/f TO i)</i + /V (Eq.8) 

T 2 = ^ hi (K m t) for j In (K m t ) > * + jV <Bq. 9) 

and: 

©-ran <E " 10) 

If the particles are spherical, K m is related to k m by: 

= (6/p*)>'""* w (Eq.ll) 

Equation 4 appears to contain five parameters, namely, a, i,j, u, and 
K m , which define the dissolution profile (W/W 0 versus time). However, 
an attempt to obtain least-squares estimates of all five parameters from 
W/ Wo versus time data may fail, because m and K m can be fused into a 
single rate parameter 

K m = e-**< m K m (Eq.12) 

which will be called the specific dissolution rate parameter. 

The uniqueness of these four parameters, a, i, j, and K %, in defining 
the dissolution profile can be seen by substituting K m = e^ m K * m into 
Eq. 4, resulting in total cancellation of p, in full agreement with the theory 
discussed previously 

W Q &k\J { Km> F0'-3 ff )-F(-i-3a) e 

(Eq. 13) 

where: 

A=sF (-i-^) forr<^re-^» (Kq.14) 



A-F [ito (O) - 2^] for )Jr > « > jjr*-**'"' 

(Eq.15) 



and: 



toxt*<gre*'' /m 



(Eq. 16) 



, The continuously recording flow-through dissolution apparatus used 
provides dissolution rate data. The fraction undissolved, W/Wq, versus 
time can be obtained by integrating these data (5). The parameters a, i, 
j, and K* m can then be estimated by nonlinear least-squares regression 
analysis using Eq. 13. However, the integrated data contain integration 
errors. The integration also tends to "smooth" the original data, so esti- 
mates of the variability of the parameters will be less reliable than if the 
original rate data were used Therefore, it is useful to derrve an expression 
for the release rate to estimate the four parameters directly from the 
original (rate) data. 

By applying Eq, 13a of Ref. 1, the following equation (similar to Eq. 
A2 in the i 



(Eq. 17) 



reference) ia obtained: 
IV Xr U3/m ~ ****** %NQn X,fl ' 0) * 



f*xmynx,,i,o)dx 



where: 



do = 

Do - e* + J° 

*, = d 0 

Ri « (K m t)^ 

R 2 'D 0 

R2=(K m tr* 



forJf m t<d 0 3/m 
forK M t>d 0 3/w 
for K M t <Do 3/m 
forff m t*D 0 3/m 



(Eq. 18) 
(Eq. 19) 
(Eq.20) 
(Eq.21) 
(Eq.22) 
(Eq.23) 



An expression for the dissolution rate of drut;, Q ™ —dW/dt, can now 
be found by differentiating Eq. 17 with respect to time (using Leibnitz's 
rule). For abbreviation, let B denote the (constant) denominator of Eq. 
17; differentiation then gives: 

-2-Q - mK m $** {x*>>» - Jft)t—»>x' 'NQnx t u,e) dx 



+ (/V' m - K m t)»Rr Win Rm) ^ (Eq. 24o) 

Before the critical tin*:* 2 , R x - d 0 and fl 2 ■ Do, so dR x /dt - dR*it - 
(» and the last two terms of the right-hand side of Eq. 24a vanish. After 
the critical time, dR*fdt is still zero, but Ri » (K m t ) m <* (according to Eq. 
21). Thus, the last term also vanishes because: 

*\* m - K m t = [(Km*)"" 3 ] 3 "" - K m t « 0 (Eq. 246) 

Equation 24o can be simplified to: 

f * (x 3 "» - K m t)<*-"x-WQn *,p,<r) dx 
JR\ 

J^°°*2N(lnw)dx 

(Eq. 25) 

under the integral sign can be expand - 



Q = mK M W 0 - 



The term (x^" - K m t)< m ~ l h 
ed: 

= x-2K3t + (Kst)*x- 1 

and a formula given previously (Eq. A3 of Ret. 1) fan be applied to express 
the integrals in Eq. 25 in terms of the function F (Eq. 5). This process 
Und* to the following expression for Q that can be evaluated more readily 



(form = 3) (Eq. 26) 
(form -2) (Eq.27) 



a The cntiuU time is the time when the first particles start disappearing in the 
dissolution prooew, £.«.. t ' (l/K w ^eV^-'*>. 
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Figure I— Illustration of a graphic method for obtaining an initial es- 
timate of the specific dissolution rate parameter, K^, from rate data 
from the dissolution of 5 mg of micronized glyburide. Equation 39 pre- 
dicts a linear relationship between VQand t for a monodisperse powder 
in which the single particles dissolve according to the cube root law (Eq. 
1, m - 3). The deviation from linearity in the graph is caused by the 
particle-size distribution effect, 

and exactly than Eq. 25: 

Q = mlf m H'o ( 1c' , ( m ; 1 )(-K m t)« » 

n-o \ n / 

F f T2-n _3n<r\ _p/ Ti-n Zna\ 
F{j - 3<r) - F(-i - 3a) 

(Eq.28) 

where Ti and T 2 are as defined previously (Eqs. 6-9). 

This expression for the dissolution rate, Q, contains five parameters 
(a, i, j, m, and K m ). However, as before (Eq. 13), only four parameters (*, 
i, j, and K m m ) are needed to define uniquely the dissolution rate pro- 
file: 

Q = mkV u W Q 'J:' ( m " l ) (-K^— » 



e «(n/m)t-i|.ift (Eq.29) 



F(j - Za) - F[-i - Ze) 
This equation is the differential form of Eq. 13. 

The quantity A is as defined previously (Eqs. 14 and 15) and Q - 0 for 
t £ (l/K m ) e$ 9/m . This expression is valuable since it allows the effective 
initial particle-size distribution parameters, <r, i, and j, to be determined 
together with the specific dissolution rate parameter, K * m , by nonlinear 
regression analysis of the dissolution rate profile (dW/dt versus time). 

It is also of interest to investigate how well the dissolution behavior 
can be described if the particle-size distribution is considered ideal, i.e., 
if i =y = «. With - 1 and F(-») = 0, Eq.29 becomes, lor an idea) 
distribution: 

Q-»^W^> < "E ,> (. m ~ , ) i-Kifih-r» 

ii'O \ n / 

X { 1 - F In (K' n t) - ^] } e*H«f»W<* (Eq. 30) 

If the distribution is considered ideal at the lower end (i ■ «) but trun- 
cated at the higher end, the expression becomes: 

«-«r.», ( J 1 ( m ; *) (-«;*)""—» 



e 91<n/m>*-l]aV2 (Eq. 31) 



FO' - 3<r) 

and if it is considered truncated at the lower end but not at the higher 
end 0 ■ »): 



[m-fi-1) 

I - A 



Discussion of the limiting case where the particle-size distribution is 
infinitely narrow* ie , o - 0, allows a better understanding of the specific 
dissolution rata parameter, K' m . It also provides a method of obtaining 
a suitable initial estimate of this parameter to use in nonlinear curve 
fitting. 

When a « i « j ? «• 0, Eq. 29 becomes: 
which can be written more simply as: 



9- 

di 



■ilRilWl-JC^^-M 



(Eq. 34) 



l-JW-3a) 



which, after integration, can be written: 

KwJ ml - K »< 

As expected, since there is no Bize distribution effect (2), this equation 
predicts that the powder will dissolve strictly according to the cube root 
or square root law. 
When the powder is monodisperse, m * In a 0 and Eq. 12 becomes: 

Km - e-M»^ °*K m - o 0 -^"K m (Eq. 36) 

Equations 11 and 36 then give, for spherical particles 1 : 

*; - a 0 - a *»<6/pir)i''"A IB - wv-M»k n (Eq. 37) 

When k m * w 0 Vm K m (from Eq. 37) is inserted in Eq. 1: 

(ST- 1 -*-' 

A comparison of Eqs. 38 and 35 shows that, for a monodisperse powder, 
the specific dissolution rate parameter, K ^ is common to both mul- 
tiparticulate (Eq. 35) and single-particle dissolutions (Eq. 38). It has the 
dimension of time -1 for both single-particle models (m • 2 and 5). 
Equation 34 can be written: 

QiHm i> „ (mlVofC.) 1 ' 0 "- 1 * - (mW 0 ) l/( " ^O™'*™-^ (ISq. 39) 

This equation can be used to obtain an initial estimate of iOrom the 
linear regression of Q^»»-*> on t. In thk way, the intercept value, 
(mWoK m )W«-"(t - 0), divided by (mWo) l ^ m ' li gives K m l ^> (Pig. 

EXPERIMENTAL 

The high precision, flow-through, continuously recording dissolution 
apparatus used for the dissolution tests was described previously (5). 
Various amounts of micronized glyburide 4 were exposed in the dissolution 
cell to pH 7.25 phosphate buffer (/ = 0.05) at a flow rate of 0.568 ml/sec 
(~ 0.149-cm/sec linear flow rate) and at a temperature of 37.5 ± 0.2°. 
Values of absorbance were read from the chart recording at 2.5-min in- 
tervals. The dissolution rate, Q, was then given as (volumetric flow rate) 
X (Bbsorbance)/(absoiption coefficient). 

The accuracy of the rate data was checked by integration. The amount 
released, calculated in this way, agreed to within *5% with the amount 
remaining in the dissolution cell at the end of the experiment (which was 
about 10% of the initial amount). The lower and upper limits for the 
nonlinear regression parameters used in the curve fitting were chosen 
a* a - 0-1, K; « 0 0,1 , i = 0-10, andy = 0 10. The initial estimate for 
the rate parameter, fC m , was obtained graphically (Eq. 41 and Fig. 2). 

RESULTS AND DISCUSSION 

For a monodisperse powder, Eq. 39 predicts a linear relationship be- 
tween VQ and t when the single particles dissolve according to the cube 
root model (m * 3) and a linear relationship between Q and t when they 
dissolve according to the square root model (m » 2). 

Signiricant deviations from linearity were observed when dissolution 
rate data for micronized glyburide were plotted in either of these ways 
(Figs. 2-6). Such deviations arise if the powder is not monodisperse or 
if the single particles do not dissolve accoidiiig fc> the sinf le-partide model 
(Eq.l), 

Under an electron microscope, the micronized glyburide used appeared 
to be considerably poly disperse Therefore, the observed deviation from 



3 The derivation fa also valid 'or aonspherical particles if the c mcept of spfomral 
approximations (8) applies. 
* Provided by Auttralian HotcUst Ltd., Sydney, Australia. 
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Figure 3— Least-squares fit of the multiparticulate dissolution model 
(Eq. 30 f va = 3) to rate data from the dissolution ofSmgof micronized 
glyburide. The model assumes that the effective initial particle-size 
distribution can be approximated by an ideal log-normal distribution 
(\ * j * a> t Fig. 1) and that the single particles dissolve according to the 
cube root model (Eq. 1, m * 3). The inset graph shows a plot of the re- 
siduals versus time in arbitrary units. 

linearity in the rate plots can be explained as a particle-size distribution 
effect if Eq. 1 is assumed to be valid. However, the rate data can also be 
explained by other single-particle dissolution models in combination with 
a size distribution effect. Conclusions about the validity of a single-par- 
ticle model can only be made when the size distribution effect is taken 
into account. 

In a previous study (5) dealing with the dissolution of 60-85-mesh 
fraction tolbutamide in relation to its particle-size distribution, the dis- 
solution of the single particles accorded well with the cube root model 
and nearly as well with the square root model. The distribution effect 
could be evaluated for tolbutamide because of the regular form of its 
particles. This was not possible for the much smaller micronized glyburide 
particles, which had very irregular shapes. 

Effective Particle-Size Distribution— Spherical particles are as- 




20 30 
MINUTES 

Figure A— Least-squares fit of the multiparticulate dissolution model 
(Eq. 29, m - 3) to rate data from the dissolution ofSmgof micronized 
glyburide. The model assumes that the effective initial particle-size 
distribution can be approximated by a log-normal distribution trun- 
cated at both ends (Fig. 1) and that the single particles dissolve ac- 
cording to the cube root model (Eq. l,m m 3) The inset graph shows a 
plot of the residuals versus time in arbitrary units. 




MINUTES 



Figure 5 — Least-squares fit of the multiparticulate dissolution model 
(Eq. 30, w 3) to rate data from the dissolution of 10 mg of micronized 
glyburide. The model assumes that the effective initial particle-size 
distribution can be approximated by an ideal log-normal distribution 
(\ m j - », Fig. I) and that the single particles dissolve according to the 
cube root model (Eq. 1, m = 3). The inset graph shows a plot of the re- 
siduals versus time in arbitrary units. 

sumed for the multiparticulate dissolution equations presented. Such 
particles are usually not encountered in pharmaceutical systems. How- 
ever, the dissolution of a nonspherical particle can, in most cases, be well 
approximated by the dissolution of a hypothetical spherical particle (8), 
particularly under conditions of isotropic dissolution. 

A previous paper showed how an equivalent spherical diameter can 
be calculated from the geometry of simple crystal forms (8). In this way, 
it w possible to obtain an effective particle-size distribution, «.*., the shape 
of the distribution of hypothetical spherical particles that approximates 
the dissolution behavior of the real particles. For tolbutamide, the ef- 
fective distribution was approximately log-normal (5). Because of the 
irregular particle shapes, it was not possible to make a priori conclusions 
about the effective particle-size distribution for glyburide. In these in- 
vestigations, the distribution was assumed to be log-normal, consistent 
with previous results (5) and the met that powders often have log-normal 
particle-size distribution (9, 10). Therefa-e, the distribution in Fig. 1 U- 




MINUTES 



Figure 6— Least squares fit of the multiparticulate dissolution kinetic 
model (Eq. 29, m - 3) to rate data from the dissolution of 10 mg of mi- 
cronized glyburide. The model assumes that the effective initial par- 
iicle-size distribution can be approximated by a log-normal distribution 
truncated at both ends (Fig. 1) and that the single particles dissolve 
according to the cube root model (Eq. 1, m - 3). The inset graph shows 
a plot of the residuals versus time in arbitrary units. 
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Tabic I— Least-Squares Estimates of Rate and Distribution Parameters Obtained from Nonlinear Regression Analysis of Data 
from the Dissolution of Micronized Glyburide, Using Various Models for Multiparticulate Dissolution Kinetics" 



mg 



Parameter 



Square Root Model (m « 2) _ 
Eq. 30* Eq. 3K Eq. 



Cube Root Model (w»3l 



Eq. 30* 



Eq. 31 c 



Eq. 29<* 



10 



> 

a 

i 
j 



min" 



0.02858 (6.39) 
0.4759 (2.93) 



0.9935 

0.02633 (7.18) 
0.4955 (3.20) 



0.9938 



0.02858 (5.12) 
0.4760 (2.70) 

5.299 
0.9935 

0.02633 (10.6) 
0.4956 (7.24) 

5.299 
0.9938 



0 02859 (6.51) 
0.4760 (3.28) 
5.136 
5.655 
0.9935 

0.02633 (9.36) 
0.4955 (6.14) 
8.912 
6.732 
0.9938 



0.02367 (4.56) 
0.6184 (2.47) 



0.9952 

0.02147 (9.45) 
0.6374 (3.29) 



0.02367 (4.44) 
0.6184 



f (4.44) 
(2.11) 



0.9957 



9.807 
0.9952 

0.02147 (4,74) 
0.8374 (2.03) 

7.060 
0.9957 



0.02449 (3.96) 
0.6315 (1.71) 
2.364 
4.172 
0.9954 

0.02176(4.03) 
0.6446 (1 89) 
4.577 
4.200 
0.9958 



a Values in parentheses arc relative standard deviations (percent). * Initial size distribution is considered ideal (i 3 Fig. 1). c Initial size 
ution is truncated at the upper end (i «=«,/< Fig. 1). * In»ial size distribution is truncated at both ends (Fig. 1). * Correlation coefficient. 

tigated. it can be concluded that the following equation; 

X {l - F [i In (K\t) - nc J J € <»«-»-Va 
which can be written in integrated form as: 



size distri- 



bution 

lustratea a log-normal approximation to the effective particle-size dis- 
tribution. The initial particle diameters, a 0 (Eq. 2), are the diameters of 
the hypothetical spherical particles that approximate the dissolution of 
the real nonspherical particles. 

If the lower and upper truncation parameters, i or are finite, the 
log-distribution is said to be truncated; otherwise, it is ideal. 

Regression Parameters— In a discussion of size distribution effects 
in multiparticulate dissolution, it was pointed out that the intrinsic 
dissolution profile does not depend on the actual size of the particles (o 0 ) 
but on the shape of their distribution (2). For the same reason, it is not 
possible to determine by regression analysis the scale parameter, m. as 
might be erroneously expected from the appearance of Eq. 2a If, in fact, 
this equation is used to obtain least-squares estimates of u and K m , these 
values would not be unique. Any other of the infinite combinations of the 
two parameters that give the same value ofK' m (Eq. 12) will, according 
to Eq. 13, result in the same fit (for the same values of i, and ;). 
Therefore, the regression analysis can only provide estimates of the di- 
mensionless distribution parameters a , i, and which define the shape 
of the initial distribution, and K" m . The scale or position of the distribu- 
tion, given by m. » bidden in the specific dissolution rate parameter, K m 
(Eq.12). 

Curve Fitting— Equation 29 was fitted by least squares to dissolution 
rate data from dissolution of 5 and 10 mg of micronized glyburide 6 . The 
fitting was done using FUNFTT, an interactive time-sharing program for 
general nonlinear regression, written by the author 6 . 

The estimates obtained for the truncation parameters, i and j (Table 
I, Eq. 29), were all larger than 2 and usually exceeded 4, indicating that 
the effective initial particle-size distribution (Fig. 1) was close to ideal. 
It was shown earner that the effect on the dissolution profile of the lower 
truncation parameter, t, is negligible (1). Simulation studies also show 
that the influence of an increase in the upper truncation parameter, 
becomes insignificant when j is larger than 2. Therefore, it is expected 
from the values of i and ; obtained using Eq. 29 that the simpler model, 
Eq. 30, which assumes an ideal distribution (i - ; - »), should fit the 
same data nearly as well. 

The values in Table I and the curves fitted in Figs. S-6 confirm this 
expectation. There does not seem to be any significant difference in the 
K^ c, or r values for the two models. The residual plots in Figs. 3-6 also 
seem to be very similar. Equation 31 (upper truncation) and Eq. 32 (lower 
truncation) also gave similar results (and, therefore, they have not been 
included). 

Choice of Mathematical Model—In agreement with the general 
principles of mathematical modeling, Eq. 30 should be considered as the 
model that best describes the dissolution of the micronized glyburide. 
It is the simplest o£ the models, containing only the two parameters K* m 
and a t and fits the dissolution data just as well as the other models that 
contain more parameters (i and j). 

It also is clear (Table I) that the multiparticulate dissolution model 
based on the cube root model (Eq. 1, m - 3) agrees best with the disso- 
lution data for either a truncated or ideal distribution using either 5 or 
10 mg of powder (Table I). Of the eight multiparticulate models inves- 



* The initial data points) « = 0 and Z5 min) was not considered in the curve 
fitting because of irregtnVilies due to initial adjustments of the solvent flow 
rate. 

« The program is available on tape from the author together with a listing of the 
subroutines specifying the multiparticulate models (Eq*. 29-31), 



(Eq.40) 



x\l-F[l]n(K%)-na\\e^-*»>» (Eq.41) 

beat describes the dissolution kinetics of the micronized glyburide. Either 
of these two equations uniquely characterizes the dissolution behavior 
in terms of the rate parameter, K* 3 , and the distribution parameter, *. 

It appears from the residual plots in Figs. 3-6 that there is a serial 
correlation between the residuals. The Durum- Watson statistic indicates 
(a < 0.05) that this correlation is significant (11). Systematic deviation 
can be caused by nonrandom experimental errors. It can also be caused 
by a significant departure from the assumed single-particle dissolution 
model (Eq. 1) or by a deviation trow log-normality. The residual values 
are, however, so small in relation to the accuracy of the experimental 
technique that the correlation seems of little importance. 

The specific dissolution rate parameter, Kl theoretically should be 
independent of the initial amount used, Wo. but this condition was only 
approximately true. Values of the rate parameter, /C obtained for W 0 
- 10 mg were consistently lower than for W 0 - 5 mg (Table I). The dis- 
solution models consider dissolution under complete sink co nd i tio ns, Le., 
conditions where there is no interaction between the dissolving particles. 
When using 10 mg of the very fine powder, it was not possible to load the 
dissolution cell with a single "layer" of particles in such a way that dis- 
solved drug from any particle did not pass over other particles. The 
slightly larger K* m values observed when leas powder was used in the cell 
agree with an expected smaller particle interaction. 

Characterization and Quality Control of Drug Powders— Previ- 
ous approaches to characterizing the dissolution properties of drug 
powders have been baaed on equations describing monodisperee systems. 
The so-called "dissolution rate constant" sometimes evaluated using such 
equations will often not charac^-.-ize the dissolution behavior because 
the size distribution effect is not en into account This ia particularly 
true for pharmaceutical systems, since these systems frequently involve 
fine, highly polydisperse powders. 

The use of nonlinear regression analysis to evaluate the specific dis- 
solution rate parameter, K* m and the distribution parameter, <r, repre- 
sents a more exact and meaningful approach. 

The properties of K* m make its interpretation particular Jv useful Ttme 
properties are best understood in relation to the concepts of time scaling 
and the intrinsic dissolution prof le (2) from which the following con- 
. .usion can be made: 

If K* m is changed by a factor a, then the tm»forccmptetedif«c4ution. 
or the time for any particular fraction to dissolve, is changed by a factor 
ofl/«\ 

TV distribution parameter, a, is a single measure of how porydisperse 



< Pwesanw*, 't it takes x nwa fora powdar to uWve, e^., 30* for a given K& 
value, then ft wffl take Aepowte^ 
of K* is doubbd (aw Pig. 1 in Ret 2). 
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a powder is, or more exactly, how much the dissolution behavior deviates 
from what would be expected if the powder were completely monodis- 
perse. A value close to zero characterizes a nearly monodisperse powder, 
while higher values indicate increasing "degrees of dispersion. Probably, 
the most important property of <r is that it is a measure of how long it 
takes the last fraction of a polydisperse powder to dissolve. For example, 
it is seen from Eq. 16 that the time for complete dissolution increases 
exponentially with For this reason, a significant correlation probably 
exists between a and systemic availability for very slightly soluble drugs 
that exhibit low systemic availability due to dissolution rate-limited 
absorption. Research in this area should be of considerable pharmaceu- 
tical interest. « 

Although the multiparticulate dissolution model (Eqs. 40 and 41) 
defining K„ and a may seem complex, the interpretation of these pa- 
rameters is simple and they can be readily obtained. The experimental 
technique used requires a high precision, flow -through dissolution ap- 
paratus that is easy to standardize [e.g., the apparatus described in Kef. 
5 or other suitable flow-through system (12)] in combination with a 
nonlinear regression program. 

The method could well become established as a routine procedure in 
quality control, and further investigation could result in improved 
standards for drug dissolution. 
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Sensitive Assay Procedure for Ethambutol 
Hydrochloride via Charge Transfer Complex Formation 
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Abstract D The charge transfer complex formation between ethambutol 
and iodine was investigated and utilized as the basis for a sensitive 
spectrophotometric procedure for ethambutol and its dosage forms. The 
solutions exhibited blue-shifted iodine bands at 293 and 360 nm. A Job's 
plot of corrected absorbance against the mole ratio of ethambutol to io- 
dine indicated a 1:2 drug-iodine ratio. At 293 nm, the absorbance was 
linear (r = 0.9998) over the 0.25-1 5-,ig/ml concentration range, but the 
concentration range for best accuracy is 1.6-5.8 Mg/ml. The method can 
be applied successfully to the analysis of commercially available 
ethambutol tablets. 

Keyphrases □ Ethambutol -spectrophotometric analysis using charge 
transfer complex formation, prepared samples and commercial tablets 
□ Spectrophotometry— analysis using charge transfer complex forma- 
tion, ethambutol in prepared samples and commercial tablets □ Charge 
transfer complex formation— spectrophotometric analysis of ethambutol 
in prepared samples and commercial tablets □ Antibacterials— etham- 
butol, spectrophotometric analysis using charge transfer complex for- 
mation, prepared samples a nd commercial tablets 

Ethambutol is a new antitubercular drug included in 
USP XIX (1). Being a relatively new drug, few procedures 
have been reported for its determination as a drug sub- 
stance or in commercial dosage forms. The official assay 
procedure involves a nonaqueous titration with perchloric 
acid titrant (1). m 

Quantitative analyses of small quantities of ethambutol 
have been based on the chelation properties of the drug 



with copper in aqueous or nonaqueous media (2, 3). A 
reineckate method and an acid-dye technique using 
bromthymol blue were also reported (4, 5). However, these 
procedures are not feasible for the analysis of micro- 
quantities of the drug in biological fluids. 

Amines and alcohols are lone-pair (n) electron donors 
and can interact strongly with sacrificial electron acceptors 
such as iodine to form charge transfer complexes (6). Re- 
cently, a study on charge transfer complexes of iodine with 
alkaloids was reported (7). This paper reports the charge 
transfer complex formation of ethambutol with iodine and 
its application to drug assay. 

EXPERIMENTAL 

Apparatus — A UV-visible double-beam spectrophotometer 1 with 
1 -cm cells and dynode voltage of 500 v (slit width 0.«M).85 mm) and an 
analytical balance 2 were used. 

Materials and Reagents— Ethambutol hydrochloride 3 and 8 X 10"* 
M iodine* (resublimed) in anhydrous chloroform 4 were used. Other re- 
agents were analytical grade. | . 

Ethambutol Base— Ethambutol hydrochloride (1 g) was dissolved 

i Beekm&a Acta V, Beckman Instruments Inc., ™}erton, Calif. 

* Metttef model H-18, MetUer toetrument* Co.. Princeton, N J. 
3 T*derle Laboratories, Pearl River, N.Y. 

* Fisher Scientific Co., Pittsburgh, Pa. 
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